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^Ch Abstract: We identify a set of "energy" functionals on the space of metrics in a given 

Kahler class on a Calabi-Yau manifold, which are bounded below and minimized uniquely 
on the Ricci-flat metric in that class. Using these functionals, we recast the problem of 
i numerically solving the Einstein equation as an optimization problem. We apply this strategy, 

using the "algebraic" metrics (metrics for which the Kahler potential is given in terms of a 
polynomial in the projective coordinates), to the Fermat quartic and to a one-parameter 
family of quintics that includes the Fermat and conifold quintics. We show that this method 
yields approximations to the Ricci-flat metric that are exponentially accurate in the degree 
of the polynomial (except at the conifold point, where the convergence is polynomial), and 
therefore orders of magnitude more accurate than the balanced metrics, previously studied as 
approximations to the Ricci-flat metric. The method is relatively fast and easy to implement. 
On the theoretical side, we also show that the functionals can be used to give a heuristic 
proof of Yau's theorem. 
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1. Introduction 

Calabi-Yau (CY) manifolds play a pivotal role in complex differential geometry, algebraic 
geometry, and string theory. Part of the reason for their importance, for both mathematics 
and physics, lies in the Ricci-flat metrics they admit. The existence of these metrics is 
guaranteed by Yau's theorem; this theorem, however, is not constructive, and, outside of flat 
metrics and orbifolds thereof, no exact solutions to the Einstein equation are known on any 
CY manifold. This state of affairs naturally leads to the question of whether it is feasible to 
compute useful numerical approximations to the Ricci-flat metrics, a question that has given 
rise to a small literature in recent years [1, 2, 3, 4, 5, 6, 7]. A few different methods have 
been put forward, all of which make crucial use of the advantages offered by Kahler geometry 
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compared to real geometry. The purpose of this paper is to introduce a new method for this 
problem, that scores well on all three of the major considerations that make such a method 
useful: ease of implementation, scale of computational resources required, and accuracy of 
the approximations obtained. The method is based on a mathematical observation that may 
be of some interest in its own right. 

The mathematical and physical motivations for solving the Einstein equation on CYs are 
discussed at length in the papers cited above, and we will not repeat them here. However, 
it is important to note that the scope of the potential applications depends strongly on the 
performance of the methods available. For example, for some problems it is sufficient to 
compute the metric at a single point in the CY moduli space, while for others it is necessary 
to scan over a part of the moduli space, requiring much faster methods. 

Any numerical method for solving a PDE must address two issues. The first is how 
to represent, or store, the function being solved for — in this case, the Kahler potential — 
for example using a real-space lattice or a spectral representation. The second is, given 
that representation, how to solve the equation, or, more precisely, how to find the best 
approximation within that representation to the exact solution. For the problem at hand, 
the principal challenges in terms of the representation are the high dimensionality of the 
manifold (usually four or six real dimensions, so that even for moderate resolutions storage 
quickly becomes a limiting factor) and its typically very complicated topology. The principal 
challenge in terms of solving the Einstein equation is its nonlinearity. 

The first methods developed for solving the Einstein equation on CYs, by Headrick and 
Wiseman [1], used a real-space lattice representation of the Kahler potential, and solved 
the equation by a Gauss-Seidel relaxation method. The method was straightforward, if not 
particularly elegant. One downside was the messiness of dealing with coordinate patches. 
The main factor limiting the accuracy of the metrics obtained was storage, due to the large 
number of lattice points required for high resolution in high dimensions (the method was 
tested on K3; storage would have severely limited the possible resolutions for threefolds). 

Subsequently, Donaldson proposed a radically different set of methods based on a spectral 
(momentum-space) representation of the metric [2]. Specifically, he advocated the use of 
so-called algebraic metrics, which had been studied from a theoretical viewpoint by Tian 
and others. We will review the definition, motivation, and properties of algebraic metrics 
in Section 3. The salient properties for our present purposes are (1) they are based on 
polynomials, and therefore very easy to work with both theoretically and computationally; 
(2) they eliminate the messiness of coordinate patches; and (3) within a certain Kahler class 
they are capable of approximating any given smooth metric — including the Ricci-flat one — 
exponentially well in k, the degree of the polynomials (k is analogous to the highest mode 
number of a finite set of Fourier modes used to approximate a function, for example, on S 1 ). 
The latter property essentially promises a form of data compression, and, when applied to 
the Ricci-flat metric, has potential to solve the storage problem mentioned above. 

1 Related methods for solving the Einstein equation on toric manifolds were discussed in [8, 9, 10]. 
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As with any spectral representation, the challenge then comes in solving the PDE. Linear 
PDEs (especially homogeneous ones) become simpler when written in momentum space, as the 
different modes decouple from each other. For nonlinear ones such as the Einstein equation, 
on the other hand, every mode is coupled to every other mode, leading to a complicated 
system of nonlinear equations. Donaldson proposed two kinds of approximate solutions to 
the Einstein equation within the framework of the algebraic metrics. The first was the so- 
called balanced metrics. For each k the balanced metric is the unique solution to a certain 
integral equation, which can be found by iterating a certain integral map. 2 Donaldson showed 
that they approach the Ricci-flat metric as k — > oo, but only as a power of k. Thus they 
do not achieve the exponential accuracy promised by the algebraic metrics. The second 
proposal, called "refined" metrics, were based on a Galerkin method: the error in solving the 
Einstein equation (more precisely, the equivalent Monge- Ampere equation) was required to 
be orthogonal to the full set of basis functions (in this case, the degree- A: polynomials). The 
refined metrics are expected to convergence exponentially to the Ricci-flat metric. However, 
efficiently computing the refined metric proved difficult. 

Now, it is well known in numerical analysis that it is far easier to minimize a well- 
behaved function than to solve a non-linear set of equations. 3 Here, by "well-behaved", 
we mean smooth, bounded below, and having no critical points other than a single global 
minimum. This is true both from a theoretical viewpoint and in terms of actual available 
algorithms; while there exist several efficient and robust numerical optimization algorithms, 
for solving non-linear equations there is essentially only the Newton-Raphson method, which 
has poor global convergence properties, and a few variants on it. So it is nearly always 
advantageous, when possible, to convert the latter type of problem into the former. In the 
case of an elliptic PDE, a well-behaved functional whose variational derivative yields the 
PDE is usually called an energy functional. Note that most elliptic PDEs do not admit 
such a variational formulation. (The utility of such a formulation is also of course clear 
from the mathematical point of view, for example for proving the existence of solutions.) 
Given an energy functional and a representation scheme (real-space, spectral, or otherwise), 
one minimizes the energy functional within the corresponding finite-dimensional space of 
functions. 4 The prototype for this method is the Rayleigh-Ritz variational method for finding 
the smallest eigenvalue of a linear operator. Especially for spectral representations, this 
strategy can vastly simplify the problem of solving the PDE. 

In order to follow this strategy for case of the Einstein equation on a CY manifold, 5 in 

Section 2 we define a family of energy functionals on the space of Kahler metrics on a given 

CY manifold, that are bounded below and minimized within each Kahler class precisely on 

2 The balanced metrics are of significant mathematical and possibly physical interest quite apart from the 
problem of solving the Einstein equation [2, 11]. 

3 For a standard discussion of this point, see e.g. section 9.6.1 of [12]. 

4 Note that the error in solving the equation will necessarily be orthogonal, in the full function space, to 
the subspace in which one is working. Hence any method of this type automatically solves some Galerkin-type 
condition. 

The possibility of following a functional-minimization strategy was briefly mentioned in [4]. 
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the Ricci-flat metric. 6 The existence of these functionals depends crucially on the magic of 
Kahler geometry and should by no means be taken for granted — comparable functionals do 
not, as far as we know, exist in the Riemannian setting (for example, the Einstein-Hilbert 
action is unbounded below, while the integrated square of the Ricci tensor has many spurious 
local minima). Note that the functionals we define are not tied to the algebraic metrics, 
and could be used to solve the Einstein equation given any spectral or other representation 
scheme. We show, incidentally, that these functionals can also be used to give a heuristic 
proof of Yau's theorem. The argument has loopholes, as we discuss, but may nonetheless be 
useful as a quick way to understand why one should expect Yau's theorem to hold. 7 

In Section 4 we combine the algebraic metrics with one of our energy functionals to define 
"optimal" metrics. We derive the Galerkin condition for the optimal metrics, and compare it 
to the one defining Donaldson's refined metrics. 

In Section 5 we report the results of applying our method to the Fermat quartic and to 
a one-real-parameter family of quintics that includes the Fermat and conifold quintics. The 
method is relatively simple to implement (requiring just a few pages of Mathematica code), 
and very efficient (running for a few minutes on a laptop computer, it yields metrics that 
are likely accurate enough for most applications). In particular, the method achieves the 
exponential accuracy promised by the algebraic metrics. As a result, the metrics found are 
far more accurate, for a given value of k, than the balanced metrics. (The high-accuracy 
metrics we obtain on the compact conifold and deformed conifold, in particular, are perhaps 
relevant for certain models of string phenomenology; see [16] and references therein). Some 
aspects of our implementation of the method are described in Section 6. 

1.1 Future directions 

The work described in this paper can be extended in many different directions. The most 
obvious is to apply the same method to CYs with less symmetry. Doing so should not require 
significant changes to the method, although attaining high-accuracy metrics will probably 
call for the use of a compiled language rather than Mathematica. As we discuss theoretically 
in subsection 3.2 and show in practice in subsection 5.2, for CYs on which the Ricci-flat 
metric contains small, highly-curved regions (e.g. near an orbifold or conifold point in the CY 
moduli space), the algebraic metrics may not provide the ideal representation of the metric. A 
real-space representation may be preferable. Nonetheless, the energy-minimization strategy 
described here should still apply. 

A second direction would be to evaluate interesting geometric quantities given the metrics 
described here, such as curvature invariants and eigenvalues of various operators, 8 and to 

6 These functionals share some of the properties of the Calabi energy [13], the Mabuchi K-energy [14], and 
its generalizations [9, 15], although they are somewhat simpler. 

7 The recent paper [15] uses a variational approach based on the Mabuchi K-energy [14] to prove a weak 
version of Yau's theorem. 

8 The papers [1, 6] discuss the calculation of scalar Laplacian eigenvalues on numerical CY metrics, with 
the second reference specializing in particular to algebraic metrics. 
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compute the metric on the moduli space using the metric on the CY. 

A third direction would be to consider the Einstein equation on non-CY Kahler manifolds. 
This could in principle be doable using the Calabi energy [13] or the Mabuchi K-energy [14]. 

A more difficult problem would be to incorporate matter fields and/or the higher-deriva- 
tive corrections to the Einstein equation that occur in string theory (as for example was done 
for the hermitian Yang-Mills equation in [3]). Here the main issue is whether the problem 
can still be expressed as the minimization of a well-behaved energy functional. Although 
we have not investigated this question, it seems likely that the answer will be yes as long 
as the geometry remains Kahler (and perhaps more generally in the context of "generalized 
geometry" ) . 

Perhaps in combination with some of the above generalizations, one could employ the 
metrics derived by our method in models of string phenomenology, as was done in [3, 5]. 

Finally, although in this paper we mainly focus on the energy functionals as calculational 
tools, it is worth exploring whether they have any direct physical interpretation, either from 
a spacetime or a world-sheet (i.e. supersymmetric sigma model) viewpoint. 



2. Energy functionals 

Let X be a compact n-dimensional Kahler manifold equipped with a fixed complex structure 
and a fixed volume form fi, which we take to be normalized, f x P- = 1- Associated to each 
Kahler form J is a volume form 

W= J T - (2-1) 
n\ 

The ratio of the two volume forms is a positive scalar that, following Donaldson [2], we denote 
rj: 

rj = ^ • (2.2) 

Yau's theorem states that each Kahler class contains a (unique) representative on which 77 is 
constant. In the following we will assume that the Kahler class is normalized, f x fij = 1, so 
the distinguished representative has 

r7=l. (2.3) 

This is a second-order elliptic PDE of -Ampere type for the Kahler potential. 

In Sections 3-6, we will assume that X is Calabi- Yau, and we will choose fi to be the 
volume form associated with its holomorphic (n, 0)-form 0: 

A=H) n fiAO. (2.4) 

The Ricci tensor for J is then R^ = —didjlnr], and the Einstein equation is equivalent to the 
Monge- Ampere equation. 

Using the metric alone, one can write down several natural and useful geometric function- 
als, such as the Einstein-Hilbert action J x \x jR and the Calabi energy J x VjR 2 [13]. (Note 



3 The recent paper [7] studied the calculation of the moduli space metric using the balanced metrics. 
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however that the Einstein-Hilbert action vanishes for any Kahler metric on a CY manifold.) 
The introduction of a background volume form, however, can be quite useful for constructing 
interesting geometric functionals, as shown by the example of Perelman's first functional [17] 
(discussed in subsection 2. 2). 10 In subsection 2.1, we will see that, using such a volume form, 
even functionals with no derivatives acting on the metric can have non-trivial properties. 
Then in subsection 2.2 we will consider a two-derivative functional similar to Perelman's. In 
Section 5 we will describe the results of numerical evaluation and minimization of both types 
of functional. 

2.1 Ultralocal functionals 

In this subsection we will show that a large class of functionals with no derivatives acting on 
the metric can yield the Monge- Ampere equation (and therefore, in the CY case, the Einstein 
equation) as their Euler-Lagrange equation. 

Let F : M + ->Sbea differentiable convex function that is bounded below, and define 

E F [J]= [ fJ,F(rj) . (2.5) 

Clearly Ep is bounded below as well. Adding a linear term to F changes Ep only by a 
constant (Ep +ar] [J] = Ep[J] + a), so without loss of generality we may assume that F(rj) 
attains its minimum at rj = 1. It follows that, in each Kahler class, Ep has a unique global 
minimum on the solution to the Monge- Ampere equation. 

We will now show that Ep has no local minima or other critical points. Under variations 
of J within its cohomology class, 

5J = iddcf), (2.6) 

the change in rj to first order is 

&-n = \rN 2 j<t>. (2.7) 

Hence the variation of Ep is 11 

5E F = f ftF'( V )5 V = ^ f lJL jF'{r l )V 2 j<j ) = l - f W V 2 /Wf (2.8) 

J X J X J X 

The Euler-Lagrange equation is thus 

V 2 J F'(r / ) = 0, (2.9) 

which, given the compactness of X and the positivity of the metric, is equivalent to the 
constancy on X of F'{rj), which in turn implies the Monge- Ampere equation. 

10 In the Kahler context one also has the Mabuchi K-energy [14], and various generalizations of it, as discussed 
for example in [9, 15]. 

n It is interesting to note that if we choose F(rf) = ijlnrj, so that Ef is (minus) the von Neumann entropy 
of the probability distribution 77 on X, then its gradient flow is the Calabi flow. 
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There is another way to think about this, which shows that the assumption of differ- 
entiability of F is superfluous. Yau's theorem establishes a one-to-one mapping between a 
given normalized Kahler class and the space of positive functions jjonl obeying j x fir] = 1. 
Given a convex function F(r/) that is minimized at 77 = 1, it is clear that the global minimum 
(and only critical point) of J x fiF(r]) on that function space — and therefore on the Kahler 
class — is rj = 1. 

Variational formulations of elliptic PDEs can often be used to establish the existence of 
solutions. Therefore, rather than assuming Yau's theorem, it is interesting to ask whether 
we can use such a functional to prove it. 12 Of course, a rigorous proof would be far outside 
the scope of this paper, but we would like to point out a simple heuristic argument, or 
"physicist's proof, using the fact that Ep is extremized precisely on solutions. The general 
principle involved is that a differentiable real function defined on an open set, that is bounded 
below and goes to infinity on the boundary of its domain, must have a minimum and therefore 
an extremum. In our case, the open set in question is the space of (smooth) Kahler metrics 
on X in a given class. The non-trivial issue is thus the behavior of Ep on the boundary of 
this space. The generic degeneration of the metric is the development of an isolated conical 
singularity, with rj going to or 00 as a power of \z\ 2 in some local coordinate system. In 
order to guarantee the divergence of Ep on such a singular metric, it is sufficient for F(rj) to 
go to infinity like (say) e n as r\ — > 00, and e 1//?7 as r\ — > 0. Of course, there are also other (non- 
generic) kinds of degeneration, in which rj goes to or 00 sufficiently slowly that such an Ep 
remains finite; indeed, for any choice of F, one can construct such a singular metric. Another 
possible degeneration is where some eigenvalues of the metric go to infinity and others go to 
zero simultaneously in such a way that 77 remains finite; then of course Ep will remain finite 
for any choice of F. 

2.2 A two-derivative functional 

Here we give an example of a functional that, like the Einstein- Hilbert action, contains two 
derivatives acting on the metric, but, unlike it, is bounded below: 

E[J] = [ pLg^d i hxr 1 djhxT 1 . (2.10) 
Jx 

A short computation shows that 

8E = a[ /V/V^d^rr 1 / 2 ) , (2.11) 
Jx 

where 4>ij = didj(j) = Sgij- Hence E is stationary if and only if r\ is constant. 
In the Calabi-Yau case, with ft = (—i) n Q A Cl, it follows immediately that 

E = - l -j^R. (2.12) 

12 The recent paper [15] uses a variational approach based on the Mabuchi K-energy [14] to prove a weak 
version of Yau's theorem. 
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This identity reveals a connection to Perelman's first functional [17], which is defined for 
a Riemannian metric but, like E, depends on a background volume form; in our present 
notation it would be written 

E P = - J ji^R + g^dahLrjdbhii^ . (2.13) 

Its gradient flow is Ricci flow supplemented by a diffeomorphism 

9ab = -Rab-VadbkLTI. (2.14) 

In the CY case, however, due to the identity (2.12), Ep vanishes identically; correspondingly, 
the gradient flow (2.14) (gtj = 0, ijij = —Vidjlnrj) is orthogonal to the space of Kahler 
metrics. The nice properties of E suggest that it may be natural to consider its gradient flow, 

m = -didj^- 1 ' 2 ) (2.15) 

(compare this to Ricci flow, = didjlnrj). 

3. Algebraic metrics 

In this section we will briefly review the definition of and motivation for the algebraic metrics. 
The discussion follows that in [2], but is hopefully translated into a more physicist- friendly 
language. More mathematically sophisticated discussions can be found in [2, 4]. 

Essentially we are searching for a nice basis of functions on our Calabi-Yau X, in which to 
expand the Kahler potential. One's first thought might be to use, for example, the eigenfunc- 
tions of some Laplacian operator. But of course we don't know those eigenfunctions, which in 
any case probably do not have any particularly simple form. Instead the strategy is to embed 
X in CF N for some N, since CP^ is a nice space that admits a highly symmetrical canonical 
metric, namely the Fubini-Study (FS) metric; as we will see, the Laplacian eigenfunctions 
with respect to the FS metric take on a very simple form. Remarkably, these eigenfunctions 
remain simple even when pulled back to X, and provide an ideal basis in which to expand 
the Kahler potential. 

We thus require X to be embeddable as an algebraic variety in CP^; furthermore, the 
Kahler class we consider on X must be the one induced from CP^ (hence it must be an 
element of the Picard lattice). 13 This is obviously a limitation; on the other hand, it includes 
many CYs of interest, both mathematically and physically. 

3.1 Algebraic metrics on CP N 

Let the homogeneous coordinates on CP^ be z a . Recall that the FS metric on CP^ has a 
Kahler potential given, in patch b (i.e. where z b ^ 0), by 



a\2 



13 The method generalizes straightforwardly to algebraic varieties in products of projective spaces, weighted 
projective spaces, etc. 
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The denominator inside the logarithm ensures that the Kahler potential transforms correctly 
on the patch overlaps. The algebraic metrics are a straightforward generalization of the FS 
metric, in which we replace the numerator with a higher-degree homogeneous polynomial. Let 



for some positive-definite hermitian matrix h. It's easy to see that the Kahler class is inde- 
pendent of the choice of k and h. The matrix elements of h are the degrees of freedom of 
the algebraic metric (except that multiplying h by a positive number just shifts the Kahler 
potential by a constant, leaving the associated metric unchanged). Hence the larger k, the 
larger the space of algebraic metrics. 

Another way to think about the algebraic metrics is in terms of the Veronese embedding of 
CP^ into CF Nk , defined by the mapping mZ. From this point of view \n(Z A h A{j Z B /\Z C \ 2 ) 
defines a generalized FS metric on CP^*, and K h is (1/k times) its pull-back to C¥ N . 

The algebraic metrics are a spectral representation in the following sense. It turns out 
that the N? functions Z A Z B / (%2 a \z a \ 2 ) k are a basis for the first k eigenspaces of the Laplacian 
(with respect to the FS metric) on CF N , in analogy to the fact that the spherical harmonics 
on S 2 can be written as polynomials in the coordinates on an embedding M 3 [2]. Therefore 
in the limit k — > oo they form a complete basis for the functions on CP^. Any metric in the 
same class as the FS metric has Kahler potential = Kf^ + <ft for some globally defined 
real function <fi. Hence when we write (3.2), we are expanding 4> (more precisely, e k ^) in a 
truncated basis of spherical harmonics, with the matrix elements of h being the coefficients. 
Note that at large k the number of degrees of freedom is N% ~ k 2N , corresponding roughly 
to k "Fourier modes" in each real direction. As usual in Fourier analysis, if eft is smooth 
then the truncated sum of spherical harmonics will approximate it exponentially well in k 
(the technical statement is that this sum will approximate it better than k u for any v [2]; in 
practice this usually means exponential convergence). 

3.2 Algebraic metrics on X 

Now let us return to X, which is sitting inside CP^ along the vanishing locus of one or 
more homogeneous polynomials in the z a . The "spherical harmonics" Z A Z B /Q^ a |z a | 2 ) fc are 
linearly independent on CP^, but not when restriced to X. However, it is very easy to see 
precisely which functions are redundant: the linear combinations of the Z A that vanish on 
X are simply the polynomials that are proportional to (any of) the defining polynomial(s) of 
X. For example, in the case N = n + 1, X is defined by a single degree N + 1 polynomial. 

14 Note that this class is not volume-normalized (does not satisfy J fj,j = 1), either on CP^ or when pulled 
back onto X. 



k be a positive integer, and Z A be the degree- A: monomials in the z a ; there are Nf. = ( t 
of them. An algebraic metric is one whose Kahler potential is 




z< 



(3.2) 
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If k > N, the space of degree k polynomials proportional to that polynomial is N k _^ N+ i^ 
dimensional. Hence the dimensionality of the quotient space is 



2 



N = 1 




N k , k<N 
N k -N k _ {N+1) , k>N 



= < 



2(k 2 + 1), N = 3 
§£:(fc 2 + 5), iV = 4 



N = 2 



(3.3) 



Note that the (real) dimensionality of the space of algebraic metrics has gone down from k 2N 
on CP^ to k 2n on X. 

The elegance and simplicity of the algebraic metrics leads immediately to two major 
practical advantages. First, one does not have to deal with coordinate patches on X, only 
those on CP^, which are quite simple. Second, they are based on polynomials, which are the 
easiest and fastest kind of function to evaluate numerically. 

A final note concerning the expected accuracy of the algebraic metric: The estimate of 
exponential accuracy mentioned above concerns the asymptotic accuracy for large k. One 
expects, however, that in order to get a good approximation, k must be sufficiently large 
that the shortest wavelength of the Fourier modes involved (roughly 1 /k) is no larger than 
the smallest characteristic size of the function being approximated. Hence we expect the 
algebraic metrics to be most useful for CYs that are not just smooth but also geometrically 
uniform, i.e. that do not have large ratios of characteristic scales. We expect this to be 
the case far from singular points in the CY moduli space. CYs that are close to orbifold 
or conifold singularities in the moduli space, for example, will contain small, highly curved 
regions. In order to accomodate the multiple length scales occurring in such cases, it may 
be more suitable to use a representation scheme whose basis functions are localized, such as 
wavelets or finite elements. 

4. Optimal metrics 

The energy functionals described in Section 2 and the algebraic metrics described in Section 3 
may be combined to yield a numerical method for solving the Einstein equation on an algebraic 
Calabi-Yau manifold X: for a given k, minimize one of the functionals Ep or E over the space 
of matrices h appearing in the algebraic metric, using an appropriate optimization algorithm. 

Out of the family of energy functionals described in Section 2, we chose the following one 
to minimize: 



(We remind the reader that here fi = (—i) n Q,AQ.) An ultralocal functional was chosen rather 
than the two-derivative functional E only because it is simpler to evaluate numerically. Within 
the family of ultralocal functionals Ep we chose F = (77— l) 2 because it is smooth, it is simple 
and therefore easy to evaluate, and it allows us to use the specialized Levenberg-Marquardt 




(4.1) 
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algorithm for minimizing sums of squares. Note that, for metrics that are sufficiently close 
to being Ricci-flat, any generic smooth F is approximately equivalent to (rj — l) 2 , since the 
distribution of rj over X is highly peaked around 1, and the first relevant term in the Taylor 
series expansion of F(rj) about rj = 1 is (rj — l) 2 (the constant and linear terms both giving 
rise to constants when integrated over X). 

We will refer to the algebraic metrics that minimize E as "optimal". From (2.8) and 
(3.2) we see that they satisfy 

For rj close to 1 this is essentially the condition that the Ricci scalar, integrated against a 
set of basis functions, vanishes, in other words it can be thought of as a Galerkin condition 
on the Einstein equation in the form R = 0. It can be contrasted with the equation defining 
Donaldson's "refined" metrics, 

which is a Galerkin condition instead on the equation rj — 1 = 0. 

As an independent check on the quality of the optimal metrics, we also evaluated the 
two-derivative functional 

E[J]= I fig Ji d l \nr,dj\nr, = -\ [ fiR (4.4) 
Jx 1 Jx 

on them. (We also experimented in a limited way with minimizing E itself.) 

Some specifics of our method for computing the optimal metrics are described in Section 
6. Although presently implemented specifically for the CYs described in the next section, it 
would be straightforward to generalize the implementation to any projectively embedded CY. 



5. Results 

We tested our method on the Fermat quartic 

4 

E(*°) 4 = 

a=l 

in CP 3 , as well as on the quintic 

5 5 

^2{z a f -5V>J]V = ° ( 5 - 2 ) 

a=l a=l 

in CP 4 , for real values of the parameter tp in the range —10 < tp < 10, which includes 
the Fermat quintic (at ■0 = 0), and the conifold (at = 1). These particular CYs were 
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chosen mainly for their simplicity and their high degree of symmetry, which vastly reduces 
the dimensionality of the parameter space, as we will discuss below. Also, this one-parameter 
family of quintics includes both singular and uniform examples, allowing us to study the effect 
of uniformity on convergence (see the discussion at the end of Section 3). Finally, the quintic 
with ip = 0.1 allowed for a direct comparison with the balanced metrics. 

5.1 Fermat quartic 




Figure 1: E versus y for the k — 2 algebraic metric on the Fermat quartic, defined by the polynomial 
(5.3), with E plotted on a normal (top) and logarithmic (bottom) scale. For comparison the FS metric 
(at y = 1 and E w 0.09) and the balanced metric (at (y,E) « (0.67,0.03); see subsection 5.2) are also 
marked. The minimum (optimal metric) is at (y,E) « (0.42,6 x 10~ 4 ). 

We will write the homogeneous degree (k,k) polynomial Z A h A gZ B appearing in the 
Kahler potential as p (so = \i\{p l / k /\z b \ 2 )). We know that the Ricci-flat metric is invariant 
under the full symmetry group of the Fermat quartic, so we may as well restrict attention to 
algebraic metrics with the same symmetry. 15 For k = 1 the only invariant polynomial is the 

15 The symmetry group is of order 3072, and is generated by permutations of the z a , multiplying each by a 
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Figure 2: E versus (x, y) for the k = 3 algebraic metric on the Fermat quartic, defined by the 
polynomial (5.4), with E plotted on a normal (top) and logarithmic (bottom) scale. For comparison 
the FS metric (at (x,y) = (0, 1) and E « 0.09) and balanced metric (at (x,y,E) (0.43,0.36,0.005)) 
are also marked. The minimum (optimal metric) is at (x,y,E) ss (0.57,0.27,2 x 10 -4 ). 



fourth root of unity (modulo multiplying all by the same root), and complex conjugating all of them. For the 
Fermat quintic, the symmetry group is analogous, and has order 150,000. For the quintic (5.2) with ip real 
but non-zero, the symmetry group is reduced: multiplying each z a by a fifth root of unity is only allowed if 
the product of the roots is 1; the order of this group is 30,000. 
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Fubini-Study one, p = J2 a \ z ' 



2 ; for k = 2 there is a one-parameter family, 



v = y f ^2\ z ' 




+ (!-») E 



a 



Z 



.o|4 . 



(5.3) 



while for = 3 there are two parameters: 




(5.4) 



Setting y = 1 in the k = 2 case, or (ar, y) = (0, 1) in the = 3 case, returns us to the FS 
metric. Figure 1 shows E versus y for k = 2, while figure 2 shows E versus (x,y) for k = 3. 
The FS metric is marked on these plots for comparison (along with the respective balanced 
metrics; see subsection 5.2). The purpose of these plots is to show that E is a well-behaved 
function, that should not present any difficulties for numerical minimization. In addition to 
being free of local minima, 16 it is well-defined and smooth for y < 0, despite the fact that the 
corresponding metrics are actually singular. (As long as there are no local minima, this is 
an advantage since the minimization algorithm may temporarily jump across to y < 0.) The 
minimum value of E for k = 2 is about 150 times smaller than the value on the FS metric, 
corresponding to a more than ten- fold decrease in the spread of rj; going from k = 2 to k = 3, 
E decreases by another factor of about 3. The minimum for k = 2 occurs at 



This k = 3 optimal metric is probably already a useful rough approximation to the Ricci-flat 
metric for some practical applications, while being computationally extremely simple. 

The number of coefficients in p is plotted in figure 3, both with and without imposing 
imposing invariance under the Fermat symmetries. The value of E for the optimal metric is 
plotted against k in figure 4 for k ranging from 1 to 17. The exponential decrease for large 
k, expected on theoretical grounds (see Section 3), is unmistakeable. Increasing A; by 1 leads 
to a decrease in E m [ n by a factor of about 8. Superimposed on this exponential decrease is 
an interesting even-odd modulation, presumably related to the even-odd modulation in the 
number of parameters evident in figure 3. In figure 5 we show the distribution of vales of r\ 
over X for the optimal metrics with k = 1, 6, 11, and 16. 

It is useful to have some independent estimate of the quality of the optimal metrics, 
other than the value of the functional E that by definition they minimize. As explained 
below equation (4.1), any generic ultralocal functional Ep is essentially equivalent to E for 

16 Although we showed in subsection 2.1 that E does not have any local minima on the full space of metrics 
in a given Kahler class, we have not proven that such minima could not appear when the functional is restriced 
to the space of algebraic metrics. 



y ps 0.42 



(5.5) 



and for k = 3 at 



(x, y) « (0.57, 0.27). 



(5.6) 
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Figure 3: Number of coefficients in the degree (k, k) homogeneous polynomial appearing in the 
algebraic metric, versus k, for the quartic. This is the number of real degrees of freedom in the 
algebraic metric (plus 1, since multiplying p by an overall constant does not change the metric). The 
top curve is without imposing any symmetry constraints; this is Nj? = 4(fc 2 + l) 2 (see equation (3.3)). 
The bottom curve is after imposing the symmetries of the Format quartic (we do not know an explicit 
formula for this number). 

metrics sufficiently close to the Ricci-flat one, so these do not serve the purpose. On the other 
hand, the functional E (4.4) is independent, since it depends on derivatives of the metric; at 
the same time it is computationally simpler to evaluate than, for example, the Calabi energy 
f /j,jR 2 . E is plotted against k for the optimal metrics in figure 6. Again the exponential 
decrease is clear, with E decreasing by a factor of about 7 for each unit increase in k. Figure 

7 shows the distribution of values of the integrand of E, I = g^di Iwqdjlnr], for k = 1, 6, 11, 
and 16. 

Another interesting possibility is to minimize E itself, rather than E. This is computa- 
tionally far more demanding, and we have not carried out a systematic investigation. However, 
our preliminary experiments indicate that the resulting metrics are not significantly different 
from those obtained by minimizing E, at least for the low values of k we have studied. Figure 

8 shows both E and E as functions of y for the k = 2 algebraic metric (5.3). The minima of 
the two functions are clearly closely aligned; in fact they are separated by less than the error 
on our numerical estimates of them. The same was true for k = 3 and 4. For k = 5 and 6 
(the highest value of k for which we successfully minimized E), the minima of E and E were 
measurably but not significantly separated. For k = 5, the value of E increased by only about 
7% going from its minimum to that of E, and similarly E increased by only about 7% going 
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Figure 4: Value of E for the optimal algebraic metric on the Fermat quartic, versus k. The exponential 
decrease is clear, accompanied by a small even-odd modulation. The line is 0.03 x 8~ k , obtained by 
fitting to the points with k — 3 through 17. 

from its minimum to that of E. For k = 6, the corresponding increases were even smaller: 
3% and 6% respectively. However, it is not clear if one can extrapolate to higher values of k 
(or other manifolds) the conclusion that minimizing E and E are effectively equivalent. 

5.2 Quintics: scanning over moduli space and comparison to balanced metrics 

The analogues of figures 3 and 4 for the quintic (5.2), with ip = —1,0,1,2,3 are shown in 
figures 9 and 10. For all cases except ij) = 1 (the conifold point), the exponential decrease 
with k is again clear, with a slightly smaller slope (in absolute value) than for the quartic. 
(The even-odd modulation seen in the quartic is also still present.) For the conifold quintic, 
however, the decrease is clearly slower than exponential. Indeed, on a log-log plot (figure 
11), the conifold points lie approximately on a straight line, with slope —4, while the Fermat 
points lie on a curve. In principle it should be possible to predict the k~ 4 convergence rate 
from the local geometric structure of the conifold singularity. 17 This is a concrete example of 
the discussion of convergence rates in subsection 3.2. These metrics are the first high-accuracy 
Ricci-flat metrics on the phenomenologically interesting [16] compact conifold and deformed 
conifold. 

Figure 12 shows how E m [ n varies with ip, for four different values of k, in the range 
—10 < ip < 10- The conifold point stands out in these plots, with a peak whose height and 

17 This turns out to be the same as the rate of convergence of the balanced metrics on a smooth CY; see 
(5.8) (note that E ~ a 2 ). 
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Figure 5: Distribution of rj — 1 for a set of 9000 points on the Fermat quartic (randomly chosen 
according to the measure fx), for the optimal metrics with k = 1, 6, 11, and 16. 



sharpness grow with k. (In particular, although the induced FS metric, k = 1, is singular, 
that singularity does not show up in E, reflecting the ultralocality of that functional — the 
induced FS metric has a curvature singularity but an everywhere finite ratio r\ of volume 
forms.) There is also a smaller and less sharp bump at ip = —1, presumably due to the 
proximity in the complex ip plane of the points ip = e ±47 ™/ 5 , which also represent the conifold 
point. (Recall that, as a coordinate on the moduli space, ip ~ e 2m ^ip.) It is also interesting 
to note that, while the Fermat quintic is a local minimium of E m [ n with respect to ip, the 
global minimum is actually past the conifold point, at ip ~ 3. However, for large (positive 
and negative) values of ip, E m i n grows rapidly, as one expects from the fact that the quintic 
is singular in the large complex structure limit. 

Figure 13 is the analogue of figure 6 for the Fermat quintic (ip = 0). Again the exponential 
decrease is clear, with a somewhat smaller slope. And figure 14 shows E plotted against tp for 
the same values of tp and k as figure 12. It shows qualitatively similar features to that figure. 
Perhaps the most interesting difference is in the k = 1 curve, representing the FS metric; 
unlike E, the two-derivative functional is sensitive to the fact that the induced FS metric is 
singular at the conifold point and at ip = oo. 
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Finally, it is interesting to compare these optimal metrics to the balanced metrics [2] . For 
comparison to the optimal metrics, the balanced metrics on the Fermat quartic for k = 2 and 
3 are marked on figures 1 and 2 respectively. The balanced metrics were thoroughly studied 
as approximations to the Ricci-flat metric in a series of papers by Douglas and collaborators 
(along with generalizations to the hermitian Yang-Mills equation and related physical issues) 
[3, 4, 5, 6]. In particular, Douglas, Karp, Lukic, and Reinbacher published enough details 
to allow a direct comparison for the quintic (5.2) with tp = 0.1 [4]. They reported that the 
following error functional, 18 

47] = ffi\r,-l\, (5.7) 

evaluated on the balanced metrics with k = 3 through 12, was given approximately by 

a » 3.1/c" 2 - 4.2/c -3 (balanced). (5.8) 

In figure 15 we plot this function, together with the value of a for the optimal algebraic metric 
for k = 1 through 10; the latter values follow, as discussed above, an exponential decrease, 
with the following fit: 

a » 0.17 x 2.2~ k (optimal). (5.9) 

18 The functional a has the same form as the ultralocal functionals Ef- However, since \r] — 1| is neither a 
smooth nor a strictly convex function, it is probably not an ideal choice for numerical minimization. 
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Figure 7: Distribution of / = g 31 di In 77 dj In 77, the integrand of 12, on a logarithmic scale, for the 
same 9000 points on the Fermat quartic and the same metrics as in figure 5. 
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Figure 8: E (top curve) and E (bottom curve) versus y for the k = 2 algebraic metric on the Fermat 
quartic, defined by the polynomial (5.3). Note that, unlike E, E goes to infinity at y = 0, where the 
metric degenerates. Their minima are separated by less than the error in their respective numerical 
estimates. 



From this plot it is clear that the optimal metrics are a significant improvement over the 
balanced ones, whether one is interested in a rough-and-ready approximation to the Ricci- 
flat metric (low k) or a high- precision one (high k). Although we have not carried out a 
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Figure 9: Analogue of figure 3 for the quintic. The top curve is without imposing any symmetries; 



the number of coefficients is 



[|fc(fc 2 + 5)] 2 . The bottom curve is after imposing the symmetries 



of the Format quintic. The middle curve is after imposing those of the quintic (5.2) with ip real but 
non-zero. 

thorough comparison of the two methods from the computational viewpoint (i.e. ease of 
implementation and computational resources required), our experience, together with that 
reported by Donaldson [2] and Douglas and collaborators [3, 4, 5, 6] suggests that the optimal 
metrics are no harder and perhaps somewhat easier to compute than the balanced metrics. 

6. Methodology 

The problem of finding the optimal metrics can be divided into three parts: 

• Constructing the relevant space of polynomials p appearing in the Kahler potential for 
a given k. 



Evaluating the functional E for a given p. 



19 



• Minimizing E over the space of p J s. 

19 While we will not describe it explicitly here, the computation of the two-derivative functional E (2.10) can 
be done following the same strategy as for E. In practice this computation is significantly more complicated. 
In particular, analytic computation of the gradient of E is prohibitively complicated, rendering its numerical 
minimization significantly more difficult. 
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Figure 10: E m \ n versus k for algebraic metrics on the quintic (5.2) for the values of ip indicated on 
the plot. The exponential decrease is clear for all values of ip except ij> = 1, which is the conifold point 
in the quintic moduli space. A linear fit of log E m i n gives E min w 0.07 x 5" fc (not plotted) for ip = 
(the Fermat quintic). 

We will explain our approach to each part in turn. All computations were performed in 
Mathematica, with a code of a few hundred lines, 20 running on an ordinary laptop. 

Some aspects of the methods described below are similar to ones employed previously in 
numerical studies of algebraic metrics [2, 3, 4, 5, 6]. 

6.1 Constructing the space of polynomials 

Our first task is to construct a basis for the relevant space of polynomials p that enter into the 
algebraic metric. This is the space of real-valued homogeneous degree (k, k) polynomials in 
the z a that are invariant under the symmetry group of X, modulo those that vanish identically 
on it. Let P be the defining polynomial of X and T its symmetry group. For the CYs studied 
in this paper, P had only real coefficients, and therefore T included as a subgroup a Z2 group 
generated by complex conjugation of the z a . (The method described here easily generalizes to 

20 This code is available for download at http://people.brandeis.edu/~headrick/physics/, along with a 
notebook illustrating its use. 
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Figure 11: Same as figure 10, except only the data for the Fermat and conifold quintics are shown, 
and k is plotted on a log scale in order to emphasize the power-law behavior of the conifold data. The 
line is the function 1.6fc -4 , obtained by fitting the conifold points to fc -4 . In contrast, the points for 
the Fermat quintic do not lie on a straight line. 



the situation where this is not the case.) Let T be the subgroup of T that acts holomorphically 
on the z a (in our case, by permuting the z a and multiplying them by (N + 1) roots of unity). 
The real-valued T-invariant (k, k) polynomials form a real vector space life. If k > N, then 
this space contains a non-trivial subspace 11^ of polynomials that vanish identically on X; 
such polynomials are of the form gP + gP, where g is a degree (k — N — l,k) T-invariant 
polynomial. We wish to construct a basis for U^/U®. 

This problem may be solved using techniques from commutative algebra, namely by 
constructing the primary and secondary invariants of T in the polynomial ring C[z a , z a \. (An 
example of such a construction is given in [5].) In practice, however, we found it easier 
to follow a less mathematically sophisticated but more direct route, and have Mathematica 
construct the necessary basis separately for each k, by: (1) constructing bases for life and 11^ 
in terms of symmetrized monomials; (2) expressing the basis elements of TI^ in terms of those 
of Ilyt; (3) using the built-in function NullSpace to construct a basis for H^/Il®. The process 
takes a few seconds. 
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Figure 12: E min versus tp for (top to bottom) k = 1,4, 7, 10 (the top curve represents the FS metric). 
(Error bars are not indicated to avoid cluttering the plot; however, the roughness of the curves, which 
is due to numerical errors, gives a good indication of their size.) 

6.2 Evaluating the energy functional 

We now turn to the evaluation, for a given p, of the functional 



This task can be subdivided into the problem of evaluating r/ at a given point, and the problem 
of integrating over X. 

First, however, let us address the small matter of the normalizations. As noted in footnote 
14, the algebraic metric is in the same class as the FS metric for any choice of k and p, but this 
class is not volume-normalized. Indeed, we do not have an analytic formula for the volume 
of X in this class, but must evaluate it numerically. Similarly, defining ft using the standard 
formula for ft for a projectively embedded CY, we do not know J x ft analytically. In practice, 
therefore, rather than calculate r\ itself we calculate the unnormalized volume ratio 
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Figure 13: Two-derivative functional E (4.4) evaluated on the optimal metrics on the Fermat quintic 
for k — 1 to 13. The line is 2.1 x 4.1 _fc , obtained by fitting to the points with k = 3 to 13. 



Defining the average with respect to £i, 



</>A - J ~Vf (6-3) 



(so that r\ = v/{v)(i), we then calculate E as the normalized standard deviation of v: 

(("-Mid 2 ). 

E " T ' (6 ' 4) 

Thanks to the simplicity of the algebraic metrics, the evaluation of v may be done "ex- 
actly" (i.e. up to round-off errors): no finite-differencing or other numerical approximations 
are necessary to evaluate the derivatives involved. Using the standard formula for Q for a 
CY embedded in CP^, we can derive the following useful formula for v. For notational con- 
creteness let us assume that we are in the patch O^+i) °f CP^, i.e. where z N+1 ^ 0, and set 
z N+1 = 1. The defining polynomial P and the mixed polynomial p are then functions of the 
N coordinates z a , where a = 1, ... ,N (and, in the case of p, their conjugates). We define 
the (N + 1) x (N + 1) matrix * o5 by 

d^p dp dp 

^ = c7^F' *^+i = ^> *™ = ^> *jv+i,JV+i=P, (6.5) 
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Figure 14: Same as figure 12, except E instead of E m i a . 
and the (N + l)-vector Q a by 

dP 

® a = ~dzP L ' Qat+i = 0. (6.6) 

We then have 

v = k 1 - N p- N det*Q0 ba Q a , (6.7) 

where ^> ba is the inverse of ^ a j. The fact that v can be computed essentially to machine 
precision is very important. Note that the distribution of values of v over X has a spread 
as small as 10 -8 in some cases (see figure 4; the spread in v is roug hly E 1 / 2 ); such precision 
would be meaningless in the presence of significant numerical errors from, for example, finite- 
differencing. 

Unfortunately, we do not know how to analytically integrate functions over X (even 
the constant function, let alone v, which is a rational function), so the integration must be 
performed numerically. We have chosen to carry it out by a Monte-Carlo method, for several 
reasons: it is simple to implement; it allows us to work directly in C¥ N , rather than having 
to construct coordinate patches on X; and it has well-understood errors. 21 However, we 
believe that it may be possible to achieve significant improvements in this area; in fact, the 
integration is the limiting step of the whole method, in terms of accuracy, efficiency, and 
memory requirements. 

21 The papers [3, 4, 5, 6] contain discussions of various Monte-Cario schemes for integration on X. 
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Figure 15: Error functional a defined in (5.7) versus k. Blue points: optimal metrics (the errors are 
negligible on this scale). Blue line: fit (5.9) to those points for k = 4 through 10. Purple curve: fit 
(5.8) reported in [4] for the balanced metrics for k = 3 through 12. 

The Monte Carlo method instructs us to pick a set of points on X randomly according 
to the measure p,, and replace the integral with an average over the points, i.e. replace ( • )jx 
with (■)«', where y! is the sum of point masses. In our implementation, the points on X 
were generated by a rejection method. We start by dividing C¥ N into patches that overlap 
only on sets of measure zero, for example Or c \ = {|z Q | 2 /|,z c | 2 < lVa 7^ c}. As in the previous 
paragraph, let us for notational concreteness work in the patch 0(jv+i) and set z N+1 = 1. 
We have the following formula for /}: 

fi = (-i) n 5 2 (P)]Jdz a f\]JdzP . (6.8) 

a p 

We then spread out the delta function a bit: 

K,c^ = H) n e(e- |P|)n^ Q A II^ ( 6 - 9 ) 

a p 

for some small parameter e (0 is the Heaviside function). The measure fi e £pN thus has 
support on a slab of thickness e about X in CP^. To choose points according to p> e £^N, we 
choose points randomly in 0(jv+i) according to the coordinate measure, and throw out those 
that do not satisfy \P\ < e. The remaining points are then projected onto X orthogonally with 
respect to the FS metric. The projected points are distributed according to the measure /t £; x 
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obtained by projecting fi e cp jv onto X; they define the measure //. This procedure introduces 
a systematic error of order e 2 , 



systematic 



(6.10) 



as well random error of order N 



r-1/2 

points' 



V- = fie,X + 0(A r p oin ts)random • 



(6.11) 



Now, a crucial fact is that these errors in approximating fx do not lead to comparable 
errors in the computed metrics. The reason is that they do not enter into the computation 
of v, which is done using the exact measure jl. Rather, they enter only in the evaluation of 
the integral, which becomes f x fj!(rj — l) 2 rather than J x fi{rj — l) 2 . The perhaps surprising 
point is that, in defining our functional E, we needn't have used ft; indeed, for any fixed 
volume form /x, the functional f x ^(rj — l) 2 is well-behaved and minimized precisely on the 
Ricci-flat metric (even more generally, so is any functional of the form J x fiF(rj), so long as 
-F'(l) = 0). Hence, strictly speaking, the errors in (6.10), (6.11) aren't errors at all; they 
just define different, equally good choices of energy functional. The main reason to integrate 
using p, is that it is canonical. While our method forces us to use a finite number of points, 
we do not wish our so-called optimal metric to depend too much on which points happen to 
be randomly selected. The criterion for choosing iVp i n ts was thus that different samples of 
points should lead to metrics that differ by a negligible amount. (The error estimates are 
discussed quantitatively in subsection 6.4.) In practice we used 22 



for all the data shown in the previous section (the larger k is, the more points are required 
to meet this criterion; a smaller number than 3000 would have sufficed for the lower values 
of k). We also chose 23 



ensuring that the e error was much smaller than the iV po i nts error. Generating a sample of 
3000 points with this value of e, on either the quartic or the quintic, required on the order 
of 10 minutes. To reiterate: to claim accuracy in the metric at the level of 10 -8 , it is not 
necessary to evaluate the integrals to this accuracy, which would have required e ~ 10~ 4 and 
-^points ~ 10 16 , which is obviously unattainable. 

22 This appears to be substantially smaller than the number of points required to accurately compute the 
balanced metrics [3, 4, 5, 6]. 

23 Except for the data required to generate figures 12 and 14, for which we chose e = 0.03, still small enough 
to ensure that the e error is much smaller than the Anoints error. There is a cost in time to choosing smaller 
values of e, since only of order e 2 points on CP^ are unrejected. For most of our experiments, the time required 
to generate the sample of points was small compared to the times required for other parts of the computations. 
When scanning the moduli space, however, that was not the case, since we required a fresh set of points at 
each point in moduli space. 



points 



= 3000 




e = 0.01, 



(6.13) 
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The minimization of E requires the repeated evaluation of v at each point for different 
choices of p. The computation of v requires the matrix ^> and the vector Q. Q is independent 
of p, and can therefore be computed once for each point and stored. Furthermore, \£ is linear 
in p, and can therefore be computed once, and stored, for each basis polynomial and each 
point. The time required to compute this data depends on k ; for a sample of 3000 points and 
the values of k studied here, it takes from seconds to hours. Storing it requires from around 
1 MB to around 100 MB. With this data in hand, evaluating E for a given p takes on the 
order of a tenth of a second. 

Because the integration is performed by evaluating the integrand at a large number of 
points in real space, the method described here is not a true spectral method, despite the 
fact that it employs a spectral representation. The amount of data associated with this 
large number of points (the data described in the previous paragraph, along with other data 
generated by the function minimization procedure) grows with k, and storage is ultimately 
the limiting factor in going to high values of k. (In principle one could compute this data, 
and even generate samples of points, on the fly, rather than storing it in RAM. However, in 
practice this is even slower than storing the data on a hard drive and swapping it in and out 
of RAM.) 

6.3 Minimizing the energy functional 

The numerical minimization of E was carried out using the Levenberg-Marquardt algorithm as 
implemented in the built-in Mathematica function FindMinimum. This algorithm is specially 
adapted to functions that are sums of squares. 24 In practice, for simplicity, the function 
minimized was actually the numerator of (6.4), as the denominator is independent of p up 
to numerical errors. The algorithm requires not only the value of v at each point but also 
its derivatives with respect to the parameters being varied (the coefficients in the polynomial 
p). A formula for these derivatives can be derived from (6.7), allowing them, like v, to be 
computed essentially to machine precision. 

In all cases the minimization was initialized on the FS metric (which is an algebraic 
metric for any k). The time required to perform the minimization ranged from a few seconds 
for low values of k to around an hour for the highest values. 

6.4 Error estimates 

The error on the minimum value of E, shown by the error bars in figure 4, was estimated by 
a bootstrap method as follows. Three independent samples of 3000 points were generated. E 
was minimized using each sample, yielding three different metrics. E was then evaluated on 
each metric using each sample, for a total of nine different values. The (logarithmic) mean 
and standard deviation of these nine values are respectively the central value and the error 
for E m i n reported on the plot. 

Finally, we give a consistency check on our calculations, taking advantage of the fact that 
the distribution of v over X for the optimal metric becomes highly peaked as k increases. By 
24 See, for example, Section 15.5 of [12] for an explanation of the Levenberg-Marquardt algorithm. 
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definition, its normalized standard deviation is E 1 / 2 , which, for the Fermat quartic, is less 
than 10~ 8 for k = 17 (the highest value for which the computation was done). The theoretical 
mean value {v)p, is independent of the metric, but the computed value (f)u' has an error on 
the order of E 1 / 2 N p( ^^ s (not counting the e-error, which is of order E l / 2 e 2 and therefore 
negligible in our computations). Figure 16 shows the normalized differences between (v)ui at 
different k and its value at k = 17, showing that (v)^ remains constant within its error, even 
as that error decreases by more than seven orders of magnitude. The error bars in the figure 
were computed using the bootstrap method described in the previous paragraph, and agree 
in order of magnitude with the parametric estimate E 1 ^ 2 N ^^ s . 
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Figure 16: Relative deviation of (v) ^ from its value at k = 17, versus k, for the optimal metrics 
on the Fermat quartic. The only difference between the plots is the scale of the vertical axis. The 
errors are estimated using a bootstrap, as described in the text. The plots show consistency, within 
the errors, in the value of (v)^. 



